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Abstract 

The continuous-time quantum walk (CTQW) on root lattice A n ( known as hexagonal 
lattice for n = 2) and honeycomb one is investigated by using spectral distribution 
method. To this aim, some association schemes are constructed from abelian group Z® n 
and two copies of finite hexagonal lattices, such that their underlying graphs tend to root 
lattice A n and honeycomb one, as the size of the underlying graphs grows to infinity. 
The CTQW on these underlying graphs is investigated by using the spectral distribution 
method and stratification of the graphs based on Terwilliger algebra, where we get the 
required results for root lattice A n and honeycomb one, from large enough underlying 
graphs. Moreover, by using the stationary phase method, the long time behavior of 
CTQW on infinite graphs is approximated with finite ones. Also it is shown that the Bose- 
Mesner algebras of our constructed association schemes (called n-variable P-polynomial) 
can be generated by n commuting generators, where raising, flat and lowering operators 
(as elements of Terwilliger algebra) are associated with each generator. A system of 
n-variable orthogonal polynomials which are special cases of generalized Gegenbauer 
polynomials is constructed, where the probability amplitudes are given by integrals over 
these polynomials or their linear combinations. Finally the suppersymmetric structure 
of finite honeycomb lattices is revealed. 

Keywords: underlying graphs of association schemes, continuous-time quan- 
tum walk, orthogonal polynomials , spectral distribution. 

PACs Index: 03.65.Ud 
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Quantum walks have recently been introduced and investigated with the hope that they may 
be useful in constructing new efficient quantum algorithms (for reviews of quantum walks, see 
[1], [2], [3]). A study of random walks on simple lattices is well known in physics(see [4]). 
Recent studies of quantum walks on more general graphs were described in [5], [6], [1], [7], 
[8]. Some of these works study the problem in the important context of algorithmic problems 
on graphs and suggest that quantum walks is a promising algorithmic technique for designing 
future quantum algorithms. 

On the other hand, the theory of association schemes [9] (the term of association scheme 
was first coined by R. C. Bose and T. Shimamoto in [10]) has its origin in the design of 
statistical experiments. The connection of association schemes to algebraic codes, strongly 
regular graphs, distance regular graphs, design theory etc., further intensified their study. A 
further step in the study of association schemes was their algebraization. This formulation was 
done by R. C. Bose and D. M. Mesner who introduced an algebra generated by the adjacency 
matrices of the association scheme, known as Bose-Mesner algebra. The other formulation was 
done by P. Terwilliger, known as the Terwilliger algebra. This algebra has been used to study 
P- and Q-polynomial schemes [11], group schemes [12, 13], and Doob schemes [14]. 

Authors in [15, 16, 17] have introduced a new method for calculating the probability am- 
plitudes of CTQW on particular graphs based on spectral distribution and algebraic com- 
binatorics structures of the graphs, where a canonical relation between the interacting Fock 
space of CTQW (i.e., Hilbert space of CTQW starting from a given site which consists of irre- 
ducible submodule of Terwilliger algebra with maximal dimension) and a system of orthogonal 
polynomials has been established which leads to the notion of quantum decomposition (QD) 
introduced in [18, 19]. In [15, 16, 19], only the particular graphs of QD type have been studied, 
where the adjacency matrices posses quantum decomposition and one can give the graph a 
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three-term recursion structure. Then, by employing the three-term recursion structure of the 
graph, one can define the Stieltjes transform of spectral distribution and obtain the corre- 
sponding spectral distribution via inverse Stieltjes transform. The QD property is inherent 
in underlying graphs of P-polynomial association schemes (for more details of P-polynomial 
association schemes, see [20], [21], [22], [11], [23]) due to the algebraic combinatorics structure 
of schemes, particularly the existence of raising, flat and lowering operators. 

Here in this work, we investigate CTQW on root lattice A n and honeycomb one by using 
spectral distribution method. In particular, we discuss the root lattice A 2 (called hexagonal 
or triangular lattice) in more details, and then generalize the results to the case of A n . To this 
purpose, first we construct some interesting association schemes from abelian group Z® n and 
finite honeycomb lattice, where in the first case, the orbits of Weyl group corresponding to the 
finite lattice, define a translation invariant (non-symmetric) association scheme on Z® n . Then, 
by symmetrization method, we construct a new symmetric association scheme, where CTQW 
is investigated on its underlying graph. In the latter case , we construct the association scheme 
from two copies of finite hexagonal lattices, where the corresponding adjacency matrix A is 
defined suitably from the adjacency matrix of finite hexagonal lattice and the other adjacency 
matrices are constructed via powers of A (in this case we have not a systematic procedure for 
construction of association scheme as in the first case). These association schemes have the 
privileges that, for large enough size of their underlying graphs, they tend to root lattice A n 
and honeycomb one, respectively. By using spectral distribution method, we study CTQW on 
these underlying graphs via their algebraic combinatorics structures such as (reference state 
dependent) Terwilliger algebras. By choosing the starting site of the walk as reference state, 
the Terwilliger algebra connected with this choice, stratifies the graph into disjoint unions 
of strata, where the amplitudes of observing the CTQW on all sites belonging to a given 
stratum are the same. This stratification is different from the one based on distance, i.e., it is 
possible that two strata with the same distance from starting site posses different probability 
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amplitudes. Then we study the CTQW on root lattice A n and honeycomb one by using the 
results of finite lattices. Moreover, by using the stationary phase method [24], the long time 
behavior of the quantum walk on infinite graphs is approximated with finite ones. In fact, 
the numerical results show that, the A 2 (honeycomb) lattice can be approximated by a finite 
hexagonal (finite honeycomb) lattice for m larger than ~ 50 (~ 60) and times t ~ 1000 
(t ~ 700). 

Another interesting property of constructed association schemes from Z® n is that, their 
corresponding Bose-Mesner algebras are generated by n commuting generators. In particular, 
the adjacency matrices are n- variable polynomials of the generators, where recursion rela- 
tions for the polynomials are given by using the structure of the association schemes. This 
property allows us to generalize the notion of P-polynomial association schemes to n-variable 
P-polynomial association schemes, where the spectral distributions associated with the gen- 
erators are functions of n variables (variables assigned to the generators). Also, we associate 
raising, lowering and flat operators with each generator via the elements of corresponding 
Terwilliger algebra. Then, by using the recursion relations associated with the Bose-Mesner 
algebra, we construct a system of n-variable orthogonal polynomials which are special cases 
of orthogonal polynomials known as generalized Gegenbauer polynomials [25], [26], where the 
probability amplitudes of the walk are given by integrals over these polynomials or their linear 
combinations. In fact, it is shown that similar to the P-polynomial case, there is a canonical 
isomorphism from the interacting Fock space of CTQW on finite root lattice A n onto the closed 
linear span of these orthogonal polynomials. Finally, we reveal the suppersymmetric structure 
of finite honeycomb lattices in the appendix. 

The organization of the paper is as follows. In section 2, we introduce briefly root lat- 
tice A n and honeycomb lattice. In section 3, we give a brief outline of association schemes, 
Bose-Mesner and Terwilliger algebras. In section 4, we give an algorithm for constructing 
some underlying graphs of so called two-variable P-polynomial association schemes and then 
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following ref. [16], we stratify the underlying graphs of constructed two- variable P-polynomial 
association schemes. In section 5, we give a brief review of spectral distribution method 
and discuss the construction of two- variable orthogonal polynomials. Section 6, is devoted to 
CTQW on hexagonal lattice and honeycomb one, by using spectral distribution method. Also, 
the asymptotic behavior of probability amplitudes of the walk at large time t, is discussed. 
Finally, we generalize the discussions of A2 to A n in section 7. The paper is ended with a brief 
conclusion together with an appendix on the suppersymmetric structure of finite honeycomb 
lattices. 

2 Root lattice A n and honeycomb lattice 
2.1 Root lattice A n 

It is well known that a Coxeter-Dynkin diagram determines a system of simple roots in the 
Euclidean space E n . The finite group W, generated by the reflections through the hyperplanes 
perpendicular to roots 



is called a Weyl group (for the theory of such groups, see [27] and [28]). An action of elements 
of the Wyle group W upon simple roots leads to a finite system of vectors, which is invariant 
with respect to W. A set of all these vectors is called a system of roots associated with a 
given Coxeter-Dynkin diagram (for a description of the correspondence between simple Lie 
algebras and Coxeter-Dynkin diagrams, see, for example, [29]). It is proven that roots of R are 
linear combinations of simple roots with integral coefficients. Moreover, there exist no roots 
which are linear combinations of simple roots ctj, i — 1, 2, n, both with positive and negative 
coefficients. The set of all linear combinations 



n(p) = p-2 




(2-1) 



n 



ciiOti I a, E Z} = ZoL.. 



(2-2) 



i=i 
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is called a root lattice corresponding to a given Coxeter-Dynkin diagram. Root system R 
which corresponds to Coxeter-Dynkin diagram of Lie algebra of the group SU(n + 1), gives 
root lattice A n . For example root system A 2 (corresponding to lie algebra of SU(3)) is shown 
in Fig.l, where the roots form a regular hexagon and a and (5 are simple roots (see Fig.l). 
This lattice sometimes is called hexagonal lattice or triangular lattice. 

It is convenient to describe root lattice, Weyl group and its orbits for the case of A n in the 
subspace of the Euclidean space E n+1 , given by the equation 

x 1 + x 2 + ... + x n+ i = 0, (2-3) 

where x±, X2,..., x n +i are the orthogonal coordinates of a point x G E n+i . The unit vectors in 
directions of these coordinates are denoted by ej, respectively. Clearly, _L ej, i ^ j. The set 
of roots is given by the vectors 

Oij = ei - ej, i^j. (2-4) 
The roots a^-, with i < j are positive and the roots 

a, = = ei - e i+1 , i — l,...,n, (2-5) 

constitute the system of simple roots. 

By means of the formula (2-1), one can find that the reflection r aij acts upon the vec- 
tor A = J2i=i m i e ii given by orthogonal coordinates, by permuting the coordinates rrii and 
rrij. Thus, W(A n ) (Weyl group corresponding to A n ) consists of all permutations of the 
orthogonal coordinates mi, m.2,..., m n+ i of a point A, that is, W(A n ) coincides with the sym- 
metric group S n+ i. The orbit 0(A), A = (m 1? m 2 , m n+1 ), consists of all different points 
(mjj,™^, ...,m in+1 ) obtained from (m 1) m 2 , ...,m„ +1 ) by permutations. 

For our purposes in this paper, we will construct an underlying graph of association scheme 
from abelian group Z® n (m > 3), such that the constructed graph can be viewed as root lattice 
A n where Z is replaced with Z m . 
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2.2 Honeycomb lattice 

The honeycomb lattice is defined by two sets of direction vectors (vectors with integer com- 
ponents), but first we should introduce the notion of odd and even vertices. A vertex is odd 
if the sum of its components is odd, otherwise it is even. The honeycomb lattice is a two 
dimensional lattice defined as follows 

Definition 4 For an even vertex, the set of direction vectors is {(1, 0), (—1, 0), (0, 1)} and for 
an odd vertex, the set of direction vectors is {(1,0), (—1, 0), (0, —1)}. 

A honeycomb structure is related to a hexagonal lattice in the following two ways 

(1) The centers of the hexagons of a honeycomb form a hexagonal lattice, with the rows 
oriented the same. 

(2) The vertices of a honeycomb, together with their centers, form a hexagonal lattice, rotated 
by the angle of n/6, and scaled by a factor l/\/3, relative to the other lattice. 

The ratio of the number of vertices and the number of hexagons is 2 (see Fig. 2). 

In section 4, we will construct an underlying graph of association scheme from two copies 
of hexagonal lattices, where the graph is equivalent to honeycomb lattice as the size of the 
graph grows to infinity. 

3 Association schemes and their Terwilliger algebra 

In this section, we give a brief review of some of the main features of symmetric association 
schemes. For further information about association schemes, the reader is referred to [9], [10], 
[11]. 

Definition 3 (Symmetric association schemes). Let V be a set of vertices, and R{ (i = 
0, ...,d) be nonempty relations on V. If the following conditions (1), (2), (3), and (4) be 
satisfied, then the pair Y = (V, {Ri}o<i<d) consisting of a vertex set V and a set of relations 
{Ri}o<i<d is called an association scheme. 
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(1) {Ri}o<i<d is a partition of V x V 

(2) i? = {(a, a) :aeV} 

(3) i2j = R\ for < i < d, where i?* = {(/3, a) : (a, /3) G 

(4) Given (a, (3) G =| {7 G V : (ck,(3) G -Ri and (7, /3) G i?,} |, where the constants 
p*j are called the intersection numbers, depend only on i,j and k and not on the choice of 
(a, /3) G R k . 

The underlying graph Y = (V,Ri) of an association scheme is an undirected connected 
graph, where the set V and i?i consist of its vertices and edges, respectively. Obviously 
replacing R\ with one of the other relations such as Ri, for i ^ 0, 1 will also give us an 
underlying graph T = (V, Ri) (not necessarily a connected graph ) with the same set of vertices 
but a new set of edges Ri. 

Let C denote the field of complex numbers. By Maty{C) we mean the set of all n x n 

matrices over C whose rows and columns are indexed by V. For each integer i (0 < i < d), let 

Ai denote the matrix in Maty(C) with (a, /3)-entry as 

f 1 if G ifc, 

(AW = («,^n (3-6) 

1 otherwise 

The matrix Ai is called an adjacency matrix of the association scheme. We then have A — I 
(by (2) above) and 

MAi = Y. V %A k , (3-7) 

k=0 

so A , Ai, A d form a basis for a commutative algebra A of Maiy(C), where A is known 
as the Bose-Mesner algebra of Y. Since the matrices Ai commute, they can be diagonalized 
simultaneously. 

We now recall the dual Bose-Mesner algebra of Y . Given a base vertex a G V, for all 
integers % define E* = E*(a) G Mat v (C) (0 < i < d) to be the diagonal matrix with ((3,(3)- 
entry 

f 1 if (a, (3) G i*, 
(£*W = (« G (3-8) 

1 otherwise 
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The matrix E* is called the i-th dual idempotent of Y with respect to a. We shall always set 
E* — for % < or % > d. From the definition, the dual idempotents satisfy the relations 

5>; = J, E*E] = b l3 El 0< i ,j<d. (3-9) 

i=0 

It follows that the matrices Eq, E{, E^ form a basis for a subalgebra A* = A* (a) of Mat v (c). 
A* is known as the dual Bose-Mesner algebra of Y with respect to a. 

Definition 4 (Terwilliger algebra) Let the scheme Y = (V, {Ri} 0<i<d ) be as in definition 1, 
pick any v £ V, and let T = T(v) denote the subalgebra of Maty{C) generated by the Bose- 
Mesner algebra A and the dual Bose-Mesner algebra A*. The algebra T is called Terwilliger 
algebra of Y with respect to v. 

Let W = C v denote the vector space over C consisting of column vectors whose coordinates 
are indexed by V and whose entries are in C. We endow W with the Hermitian inner product 
(, ) which satisfies (u,v) = u l v for all u, v e W , where t denotes the transpose and - denotes 
the complex conjugation. For all f3 e V, let \f3) denote the element of W with a 1 in the (3 
coordinate and in all other coordinates. We observe : (5 e V} is an orthonormal basis 
for W. Using (3-8) we have 

Wi = E*W = span{\(3) : (3 G V, (a, (3) e Rt}, 0<i<d. (3-10) 

Now using the relations (3-9), one can show that the operator E* projects W onto Wi, thus 
we have 

w = w e Wi © • • • e w d . (3-11) 

In [16], CTQW on some special kinds of underlying graphs of P-polynomial association 
schemes has been investigated. It is shown in [20] that in the case of P-polynomial association 
schemes, A4 = pi(A)(0 < i < d), where pi is a polynomial of degree i with real coefficients. 
In particular, A generates the Bose-Mesner algebra. Moreover, for a P-polynomial scheme, 
there is a quantum decomposition for adjacency matrix of the underlying graph, where in 
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[16], this property has been employed for investigation of CTQW via spectral distribution 
associated with adjacency matrix. In fact, for P-polynomial schemes a quantum decomposition 
for adjacency matrix can be defined by the following lemma 

Lemma (Terwilliger [11]). Let T denote an underlying graph of a P-polynomial association 
scheme with diameter d. Fix any vertex a of T, and write E* = E*(a) (0 < i < d), A x = A 
and T= T (a). Define A~ = A~(a), A = A (a), A+ = A + (a) by 

A' = E E*_,AEl A = £ E*AE*, A + = £ E* +1 AE*. (3-12) 

i=l i=l i=l 

Then 

A = A + + A~ + A , (3-13) 
where, this is quantum decomposition of adjacency matrix A such that, 

(A-) 1 = A + , {A ) 1 = A°, (3-14) 

which can be verified easily. 

Note that the above lemma is true only in the cases of P-polynomial association schemes. 
In this paper we will construct some underlying graphs of association schemes for which the 
corresponding Bose-Mesner algebras are generated by n commuting operators. Hereafter, 
we will refer to these types of association schemes as n-variable P-polynomial association 
schemes. As a generalization of the above lemma to n-variable P-polynomial association 
schemes, one can define raising, lowering and flat operators as in (3-12) with respect to each 
generator of Bose-Mesner algebra. In particular, for association scheme derived from Z m x Z m , 
the corresponding Bose-Mesner algebra is generated by two commuting operators A z and A 2 
(Az = A\), i.e., A k i = Pm(A z , A 2 ), where p ki is a polynomial of degree k+l with real coefficients. 
The raising, lowering and flat operators are defined as in (3-12) with respect to each generator 
of Bose-Mesner algebra. Explicitly we have 

At ■■= E e; +1 a z e* At := E e; +1 a z e;, 

i i 
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Aj:=(At)< A;:=(AtY and 

A\ := £ E*A z E* A\ : = £ E*A- Z E*. (3-15) 

i i 

Similar to P-polynomial association schemes, we have 

A z = At + A; + A° z , A- z = Aj + A~ z + A% (3-16) 



4 construction of some translation invariant association 
schemes 

In this section, we construct two types of finite underlying graphs of association schemes from 
finite abelian group Z m x Z m (m > 3) and two copies of finite hexagonal lattices, such that in 
the limit of the large size of the graphs, the underlying graphs tend to infinite graphs on root 
lattice A 2 and honeycomb one, respectively. We will show that the corresponding Bose-Mesner 
algebras are generated by two commuting operators, in particular all elements of Bose-Mesner 
algebras are two-variable polynomials of the generators. We will refer to these schemes as 
two- variable P-polynomial association schemes. To our purpose, first we give some definitions. 

Definition 6 Let A be a finite multiplicative abelian group and R = {Rq, ■ R r } a collec- 
tion of r + 1 distinct relations on A forming a partition of the cartesian power A 2 . If (x, y) G Ri 
implies (ax, ay) G Ri for all a G A and i — 0,1, r, then P is called translation invariant. 

Definition 7 A partition P = {P , P r } of an abelian group A is called a blueprint [9] if 

(1) P = {e} (e is the identity of the group), 

(2) for % = 1, ...,r, if x G Pi then x~ l G P (i.e., P = P' 1 ), 

(3) there are integers such that if y G P& then there are precisely elements x G P such 
that x~ l y G Pj. 

Now let A be an abelian group, and P = {P ,...,P r } be a blueprint of A. Let T(P) = 
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{Rq, R r } be the set of relations 

R l = {(x,y)eA 2 | x^yeP*}, (4-17) 

on A. One can notice that, if Pj is a generating set for the group A, then the underlying 
graph T = (A,Ri) is called a Cayley graph on A. From (4-17), it can be easily seen that 
R = {R , ...,R r } forms a translation invariant partition of A 2 , where Rq is diagonal relation. 
Also from condition (2) in definition 7, (x,y) G Ri implies that (y,x) G Ri, i.e., R~ l = Ri. 

4.1 construction of two- variable P-polynomial association schemes 
from Z m x Z m (m > 3) 

First we choose the ordering of elements of Z m follows 

V = {e, a, a m ~\ b, ab, a m_1 6, 6 m_1 , a6 m_1 , o^ft™- 1 }, (4-18) 

where a m = b m = e. We use the notation (k, I) for the element a k b l of the group. Clearly, 
(k, l)(k', I') = (k + k',l + V) and (k, l)' 1 = (-k, -I). Then the vertex set V of the graph will 
be {(k, I) : k,l G {0, 1, m — 1}}. Now we choose generating set 

P 10 = {(1, 0), (0, 1), (m - 1, m - 1)}, (4-19) 

for Z m x Z m . With this choice, we obtain Cayley graph T = (V, Rio), where V = Z m x Z m and 
Pio is defined by (4-17). Now, we obtain the orbits of Weyl group S3 (all possible permutations 
of (1, 0), (0, 1) and (m — 1, m — 1)). Then, the orbits 

P kl :=0((k,-l)), (4-20) 

form a partition P for Z m x Z m , where P o = {(0, 0)} (in this case, P is called homogeneous). 
Therefore, by using (4-17), we obtain a coloring for the Cayley graph (V, Rio) (with Pi 7^ Rw )- 
Clearly, for the relations R k i defined by (4-17) we have, irRkin^ 1 = R M for every n G ^3, i.e., 
((xi,x 2 ), (1/1,1/2)) G R k i iff (ir(x 1 ,x 2 ),ir(y 1 ,y 2 )) G R k i- 



continuous-time Quantum walk 14 

Moreover, since any product of two orbits Pk x k 2 an d Pi x i 2 is invariant under symmetric group 
S3, the set of orbits (consequently the set of relations Rk±k 2 ) is closed under multiplication. 
Also, if we use the notation i = (11,12), j = (31,32) and k = (ki,k 2 ), it can be easily shown 
that, for ((x,x'), (y,y')) € Rkik 2 , the intersection number 

Pij = I {(*,*') : ((x,x'),(z,z')) e R ili2 , ((z,z'),(y,y')) e R jlh }\ 

= \{(z,z') : (z-x,z'-x')eP ili2 , (y-z,y'-z')eP nn }\, (4-21) 

is independent of the choice of ((x,x'),(y,y')) G R klk2 - Therefore, the relations R k i define 
an abelian association scheme (not necessarily symmetric) on Z m x Z m , where in the regular 
representation of the group, for the corresponding adjacency matrices we have 

A k ,i = E 9- (4-22) 
geP k ,i 

From (4-20) and (4-22), it follows that the adjacency matrices satisfy the following recursion 
relations 

AioA kt i = A k+ i t i + A k j_i + Ak-1,1+1 , 

AoiA k j — A k -ij + A k j + i + A k+ i jt _i , (4-23) 

where, A 0Q = I, A w and A m = A\ are the first adjacency matrices. In fact, the following two 
matrices 

A z := S! + S 2 + (S 1 S 2 )- 1 and A- z := (A z f, (4-24) 

generate the whole Bose-Mesner algebra of above constructed association schemes. In partic- 
ular, Aki = Pki(A z , A z ), where pki is a polynomial of degree k + 1 with real coefficients. We will 
refer to these types of association schemes as two- variable P-polynomial association schemes. 

We illustrate the construction of underlying graph in simplest case m = 3 in the following 
Example: case m = 3 

From (4-20), the orbits of Weyl group S3 are obtained as 

Poo = {(0,0)}, P 10 = O((l,0)) = {(l,0),(0,l),(m- l,m- 1)}, 
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Poi := O((0, -1)) = {(2, 0), (0, 2), (1, 1)}, P u := 0((1, -1)) = {(1, 2), (2, 1)}. (4-25) 

Now by using (4-17), one can obtain the relations i4ifc 2 > f° r (ki, k 2 ) G {(0, 0), (1, 0), (0, 1), (1, 1)}. 
Also, it can be verified that T(P) = {R^k?} is an abelian association scheme. The basis of 
Bose-Mesner algebra and dual Bose-Mesner algebra are 

A 00 = J 9 , A w = S 1 + S 2 + (S^) 2 , A 01 = Sf + S 2 2 + 5x5a, A u = + S 2 S 2 and 

E* = E' *®E' *, E* w = E' *®E[*+E' 1 *®E , *+E' 2 *®E , 2 *, E* Q1 = E' *®E' 2 *+E 2 *®E' *+E' 1 *®E' 1 * , 

E n * = E[* ® E' 2 * + E' 2 * ® E[*, (4-26) 

respectively, where 

(E'*) yy = 6 yi , 2 = 0,1,2. (4-27) 
The adjacency matrices are written in terms of A z and A 2 as follows 

A 00 = I, A W = A Z , A 01 = A- Z , A u = ^{A z A z -3). (4-28) 

One can notice that, the set {^i, S 2 } is a generating set for Z m x Z m , i.e., the elements of 
the group in the regular representation are of the form (k, I) = SiS l 2 , for k, I G {0, 1, m — 1}. 
If we represent Sj as Sj = e 2mx ^ m , Xj G {0, 1, ...,m — 1}, then we have SjS k = e 2m {x 3 +x k )/ m ^ 
so the multiplication in the generating set {Si,i = 1,2} is equivalent to the addition in the set 
{xi,i = 1,2}. In the additive notation, A z is written as 

— e ^ixi/m _|_ e 2Tvix 2 /m _|_ e -2wi(x 1 +x 2 ) /m (4-29) 

so, clearly {(xi,x 2 ,X3 = —{x\ + x 2 )) : Xi G Z m } is a finite sequence of triples such that in 
the limit of large m tends to the root lattice A 2 . 

4.1.1 finite hexagonal lattice 

The underlying graphs of two-variable P-polynomial association schemes constructed in pre- 
vious section are directed graphs since the relation R w is non-symmetric. In this section, in 
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order to obtain undirected (symmetric) underlying graphs of two-variable P-polynomial asso- 
ciation schemes, we symmetrize the above constructed graphs of previous section. To do so, 
we choose a suitable union of the orbits such that the new partition Q is symmetric in the 
sense that Qu = Q^i, for all (k, I). In another words, we construct a blueprint from partition 
P, by symmetrization. Such a symmetrization conserves the property of being association 
scheme, because the union of the orbits is still invariant under the action of symmetric group. 
In appendix A of [16], such a symmetrization method is used for group association schemes. 

Therefore, we construct the new underlying graph of association scheme, by choosing the 
generating set Qio as follows 

Q 10 = P 10 U P 01 = {(1, 0), (0, 1), (1, 1), (m - 1, 0), (0, m - 1), (m - 1, m - 1)}. (4-30) 

With this choice, the adjacency matrix of underlying graph is 

A = A Z + A- Z , (4-31) 

where, A z and A? are defined in (4-24). Clearly, the new graph can be viewed as finite 

hexagonal lattice. In the following, we give the symmetric partition Q and corresponding 

adjacency matrices of underlying graph for m = 3. 

Example: case m = 3 

Using (4-25), the new partition Q is given as 

Qoo = {(0, 0)}, Q w = {(1, 0), (0, 1), (1, 1), (m -1,0), (0, m - 1), (m - 1, m - 1)}, 

Qn = Pi,i = {(l,2),(2,l)}, (4-32) 

and for the adjacency matrices, we have 

A 00 = I, A w = S ± + S 2 + (S 1 S 2 ) 2 + (S 1 ) 2 + (S 2 ) 2 + S 1 S 2 , A u = S 1 S 2 2 + SfS 2 . (4-33) 

Clearly the new constructed graphs are also underlying graphs of two- variable P-polynomial 
association schemes. For example in the case of m = 3 we can write 

Aoo = l, A W = A Z + A- Z , A u = ^(A z A- z -3), (4-34) 
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In section 6, we will investigate the behavior of CTQW on these undirected graphs via 
spectral method, so we need to know the spectrum of adjacency matrix A. The spectrum of 
A z in (4-24) can be easily determined as 

Zij = J + u j + uj~ {i+j) , oo = e 2ni/m ; i, j e {0, 1, m - 1}. (4-35) 

Then, from (4-35) and that the spectrum of A z is complex conjugate of the spectrum of A z , 
one can calculate the spectrum of A as follows 

Am = Zki + z* ki = 2(cos(27r/c/m) + cos(27r//m) + cos(2-7r(/c + l)/m)). (4-36) 

4.2 construction of association scheme from two copies of hexagonal 
lattice 

We extend the group Z m x Z m by direct product with Z 2 and obtain Z 2 
vertex set for underlying graph of association scheme that we want to construct. As regards 
the argument of section 2, we know that, finite honeycomb lattice is equivalent to two copies of 
finite hexagonal lattice (see Fig. 2), therefore we define the adjacency matrix A corresponding 
to finite honeycomb lattice, such that A 2 gives us Ahexagonai, the adjacency matrix of finite 
hexagonal lattice. That is we have 

A = a + ® B* + <t_ ® B, (4-37) 
where, B = I + S l + S 2 \ Clearly B l B = BB l = Si + S 2 + S X S 2 + Sf 1 + S 2 X + (Si^)- 1 = 

Ahexagonai ■ 

By computing the powers of adjacency matrix A, one can construct other adjacency ma- 
trices associated with an association scheme (not necessarily P-polynomial). Unfortunately, 
in this case we are not able to construct the association scheme via a systematic procedure 
based on group theoretical approach as in the case of finite hexagonal lattice (this shows the 
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preference of group theoretical approach). Also, it should be noted that, in this case the asso- 
ciation scheme is defined in terms of matrices (see third definition of an association scheme in 
[9]). For example we give the adjacency matrices of Bose-Mesner algebra for m — 3. 
case m = 3 

The adjacency matrices of Bose-Mesner algebra are written as 



A Q = I 2 ®I 9 , Ai = cr + ® B l + cr_ ® B, A 2 = I 2 ®A trL , 
A 3 = a + (g)(S 1 +S 2 2 +S 1 S 2 +(S 1 S 2 ) 2 +SlS 2 +S 1 S 2 2 )+a^(g)(Sf+S 2 + (S 1 S 2 ) 2 +S 1 S 2 +S 2 1 S 2 +S 1 S 2 2 ), 



One can see that Ai for i — 1, ...,4 are symmetric and X^=o^ = ■As- Also it can be verified 
that, {Ai, i = 1,...,4} is closed under multiplication and therefore, the set of matrices 
A , ...,A± form a symmetric association scheme. We give only the following multiplications of 
adjacency matrices, where we will use them later 



We will denote the graph constructed as above by r s . It is notable that, in the limit of large m, 
the graph T s can be viewed as a graph with vertices belonging to honeycomb lattice. In fact, 
starting from site e of a hexagonal lattice, the generators cl, cSf 1 an d cS 2 (c 2 = 1), generate 
the honeycomb lattice (See Fig. 2). From Fig.2 one can see that, moving on the honeycomb 
lattice by steps of length two, is equivalent to moving on hexagonal lattice by steps of length 
one. 

One can notice that, the graph T s is a bipartite graph and has supersymmetric structure 
in the sense of Ref. [30], where we discuss the supersymmetric structure of T s in appendix. 



A 4 = I 2 ®(S 1 Sl + S*S 2 ). 



(4-38) 



A\ = 3A + A 2 , A t A 2 = 2A X + 2A 3 , A ± A 3 = 2A 2 + 3A 4 , A X A± = A 3 . 



(4-39) 
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4.3 Stratification 

In this section, first we recall some of the main features of stratification for underlying graphs of 
association schemes (see for example [16]) and then stratify the underlying graphs of association 
schemes constructed in previous subsections. 

Let V be the vertex set of an underlying graph T of association scheme. For a given vertex 
a E V, the set of vertices having relation Ri with a is denoted by Ti(a) = {(3 E V : (a, (3) E 
Ri}. Therefore, the vertex set V can be written as disjoint union of Ti(a) for % = 0, 1, 2, d 
(where, d is diameter of the corresponding association scheme), i.e., 

V=(jT t (a). (4-40) 

i=0 

Now, we fix a point o E V as an origin of the underlying graph, called reference vertex. Then, 
the relation (4-40) stratifies the graph into a disjoint union of strata (associate classes) Ti(o). 

With each stratum Tj(o) we associate a unit vector in l 2 (V) (called unit vector of i-th 
stratum) defined by 

l&> = 4= E \a)£E*W, (4-41) 
V a * aer x (o) 

where, |a) denotes the eigenket of a-th vertex at the associate class Tj(o) and Oj = |Tj(o)|. 
For < % < d the unit vectors \<pi) of Eq.(4-41) form a basis for irreducible submodule of 
corresponding Terwilliger algebra with maximal dimension denoted by W ([11], Lemma 3.6). 
The closed subspace of 1 2 {V) spanned by is denoted by A(G). Since becomes a 

complete orthonormal basis of A(G), we often write 

A(G)=E©C|0 4 >. (4-42) 

i 

In the graphs constructed from Z m x Z m , the vertex set is V — {(k, I) : k, I E {0, 1, m — 
1}}. Therefore, for a given vertex (to, n) E V, T k i((m,n)) = {(m',n') : ((to, n), (m',n')) E Rki} 
is equivalent to 

r kl ((m, n)) = {(to', n') : (m' - to, n' - n) E 0((k, -/))}, (4-43) 
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where, 0((k, —I)) denote the orbits of Weyl group corresponding to the finite lattice. Now, we 
fix the vertex (0, 0) G V as an origin of the underlying graph, called reference vertex. Then, the 
relation (4-40) stratifies the graph into a disjoint union of associate classes r H ((0, 0)). Then, 
the unit vectors (4-41) are written as 

\<f>ki) = -^= E K n >> ( 4 ~ 44 ) 



akl (m,n)er fci ((0,0)) 

where, a kt = |r H ((0,0))|. In section 6, we will deal with the CTQW on the constructed 
underlying graphs, where the strata {\<pki)} span a closed subspace (irreducible submodule 
of corresponding Terwilliger algebra with maximal dimension called walk space), where the 
quantum walk remains on it forever. 

For reference state |0 O o) = 1 00) we have 

A k i\<foo)= E K n >- ( 4 ~ 45 ) 

(m,n)er w ((o,o)) 

Then by using unit vectors (4-44) and (4-45) one can see that 

A M \<f) 00 ) = Vowl^fci)- (4-46) 



In the case of finite honeycomb lattice, we have two sets of odd and even vertices, i.e., 
V = V Q +V e , where V Q is the set of odd vertices defined by V Q = {(1; k,l) : k,l G {0, 1, m— 1}} 
and V e is the set of even vertices defined by V e = {(0; k,l) : k,l e {0, 1, m — 1}}. We define 
stratum r,(«; k, I) as 

Ti(u; k, I) = {(v, k', I') : (A)(( u; M)>;fc',*'))=i}> (4-47) 

where, u, v G 0, 1 and k, I, k', I' G {0, 1, m — 1}. Now, we fix the vertex (0; 0,0) 6 V as an 
origin of the underlying graph. Then, the relation (4-40) stratifies the graph into a disjoint 
union of associate classes Fj((0;0,0)) and the relations (4-44), (4-45) and (4-46) are satisfied 
by replacing 1^(0, 0) with ^((0; 0, 0)). 
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One should notice that, these types of stratifications are different from the one based on 
distance, i.e., it is possible that two strata with the same distance from starting site posses 
different probability amplitudes. 

5 Spectral distribution 

In this section we give a brief review of spectral distributions for operators. Although the 
spectrum of underlying graphs on which we study CTQW, is easily evaluated and so CTQW 
can be investigated without spectral distribution approach, but in the limit of large size of the 
finite graphs, the best approach for calculating expected values of adjacency matrices is spectral 
distribution one. As we will see later, based on spectral distribution, one can approximate the 
behavior of the CTQW on infinite graphs with finite ones via stationary phase approximation 
method. Also, the spectral distribution approach is the best method for studying central limit 
theorems for quantum walks on graphs, see for example [31], [32]. 

In [16] and [17], CTQW on underlying graphs of QD type is investigated via spectral 
distribution, where the spectral measures associated with the adjacency matrices are single 
variable. In the case of n-variable P-polynomial association schemes, spectral measures are 
n-variable functions. Therefore, in the following we generalize the discussions in [16] and [17] 
to the case of n-variable P-polynomial association schemes. 

It is well known that, for every set of commuting operators (A Z1 , A Zn ) and a reference 
state \<f>o), it can be assigned a distribution measure u as follows 



where E(z 1 , z n ) = Y^i\ u i"' ){ u i"' \ is the operator of projection onto the common 
eigenspace of A Zl ,...,A Zn corresponding to eigenvalues zi,...,z n , respectively. Then, for any 
n-variable polynomial P(A ZI , ...,A Zn ) we have 




(5-48) 




(5-49) 
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where for discrete spectrum the above integrals are replaced by summation. Using the relations 
(5-48) and (5-49), we have 

((f) \P(A zl , ...,A Zn )\<f) ) = J ... J P(z 1 , z n )ii(z 1 , z n )dzi...dz n . (5-50) 

The existence of a spectral distribution satisfying (5-50) is a consequence of Hamburgers the- 
orem, see e.g., Shohat and Tamarkin [[33], Theorem 1.2]. 

Actually the spectral analysis of operators is an important issue in quantum mechan- 
ics, operator theory and mathematical physics [34, 35]. As an example fj,(dx) = \ip(x)\ 2 dx 
(fi(dp) = \ip{p)\ 2 dp) is a spectral distribution which is assigned to the position (momentum) 
operator X(P). Moreover, in general quasi-distributions are the assigned spectral distribu- 
tions of two hermitian non-commuting operators with a prescribed ordering. For example the 
Wigner distribution in phase space is the assigned spectral distribution for two non-commuting 
operators X (shift operator) and P (momentum operator) with Wyle-ordering among them 
[36, 37]. 

5.1 construction of orthogonal polynomials 

As regards the arguments of section 4, the Bose-Mesner algebra corresponding to two- variable 
P-polynomial association scheme derived from orbits of Wyel group corresponding to finite 
hexagonal lattice, is generated by A z and A z defined by (4-24). We assign the variables z and 
z to A z and A z , respectively. Then, in the limit of the large size of the underlying graph, the 
recursion relations (4-23) define a set of two- variable polynomials pk,i with the first polynomials 
and recursion relations as follows 

-Po,o — 1) Pi,o — z , Pq,i — z ? 
zPk,i = Pk+i,i + Pk,l-1 + Pk-l,l+l i 

zPk,i = Pk,i+i + Pk-i,i + Pk+i,i-i ? (5-51) 
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where, the polynomials P m , n in (5-51) are orthogonal with respect to the constant measure 
fx(xi,x 2 ) = 1 (where, z = e lXl + e" 2 + e ~ t{ - Xl+X2 ">) , i.e., we have 

Pm,nPm' ,n'dXid,X2 = S m ,m'8n,n'- (5-52) 

From (4-23) and (4-46), it can be seen that, there is a canonical isomorphism from the 
interacting Fock space of CTQW (irreducible submodule of Terwilliger algebra with high- 
est dimension) on the symmetric underlying graphs of two- variable P-polynomial association 
schemes derived from Z m x Z m (finite hexagonal lattice) onto the closed linear span of the 
orthogonal polynomials generated by recursion relations (5-51). In fact, the adjacency matri- 
ces of non-symmetric association schemes constructed from Z m x Z m in section 4, are equal 
to polynomials P mjTl (A z , A 2 ) and the symmetrization of the association schemes is equivalent 
to realification of two-variable polynomials P mjn (z,z). Therefore, the adjacency matrices of 
symmetric association schemes derived from Z m x Z m , are of the form P m , n (z, if P m , n (z, z) 
is real or of the form P m , n (z, z) + P mtn (z, z) if P m , n (z, z) is complex. 

It should be noted that, in the case of finite hexagonal lattice, the polynomials P k ,\ are 
not independent. Also, it can be shown that, these polynomials can be derived by using the 
raising operators and Af defined by (3-15) corresponding to symmetric underlying graphs. 
In the following, we list the strata and corresponding polynomials in the order of their first 
appearances as 

|0 O ) — P ,0 

|0 1)O ) = A+|0 O ) - P,o, |0o,i) = At\<M - Po,i, 

|0 2 , O ) = (A+) 2 |0 O ) -> P 2>0 , 10!,!) = AtAt\<f>o) - Pl,l, 100,2) = (^) 2 |0 O ) - P ,2, 

|03, o ) = (A 2 + ) 3 |0 O ) - P s ,o, |0 2 ,i) = {AtfAi^) -> P 2 ,i, |0 li2 ) = (A+) 2 A+|0 o ) -> P 1)2 , |0o, 3 ) = (^) 3 |0o) 

(5-53) 

For the sake of clarity, we construct the polynomials in the simplest case m — 3. 
Example: case m = 3 
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In this case, we have A+ = £i=o E"* +1 A Z E"* and A± = £i=o E"* i+1 A Z E"*, where A z and A- z 
are given by (4-24) and the basis of dual Bose-Mesner algebra is given by 

E"* = E*, E"\ = El + El E"* 2 = El (5-54) 

where, E* for i = 0,1,2,3 are given in (4-26). Now, by using and Af, we obtain the 
following states 

\(f> lfi ) = A+|00) = |02) + |20) + |11) = A g \00), 
|0 O>1 ) = ^|00) = |01> + |10) + |22) = A- z \00), 
= AtAt\00) = 3(|12) + |21)) = {A Z A- Z - 3/)|00). (5-55) 
Therefore, the two-variable orthogonal polynomials associated with 101,0), |0o,i) an d |0i,i) are 

Pi,o = z, P ,i = z and P hl = zz - 3, (5-56) 

respectively. Moreover, the polynomials P m , n are special cases of orthogonal polynomials known 
as generalized Gegenbauer polynomials [25], [26]. These polynomials also can be derived from 
solving the schrodinger equation for special case of completely integrable quantum Calogero- 
Sutherland model of A n type with constant potential, which describes the mutual interaction of 
N = n+1 particles moving on the circle. The coordinates of these particles are Xj, j — 1, N 
and the Schrodinger equation reads as 



N Q2 

H* = EV, P = --A, A = J2 g—j. (5-57) 



The ground-state energy and (non-normalized) wavefunction are 

P = 0, * (xi) = 1. (5-58) 
The excited states depend on an n-tuple of quantum numbers m = (mi, 777.2, ■■■,m n )' 

H* m (x i )=E m * m ,, P m = 2(A,A), (5-59) 
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where A is the highest weight of the representation of A n labeled by m, i.e., A = Y^l=\ m i e i 
and ti are the fundamental weights of A n . In fact, the eigenfunctions ^ m are solutions to the 
Laplace equation 

-A^ m = E m m m . (5-60) 
Let us restrict ourselves to the case A 2 . If we change the variables as 

Zl = e 2ixi + e 2iX2 + e 2ix \ z 2 = e 2i(xi+X2) + e 2i{x2+xs) + e 2i{x3+Xl \ z 3 = e 2i ^ xl+X2+X3 \ (5-61) 

then, in the center-of-mass frame (J2i x i — 0), the wavefunctions depend only on two variables 
chosen as z = z\ and z = z 2 (in this case, z% = 1). With this change of variables and 
using normalization for \l/ m such that the coefficient at the highest monomial is equal to one, 
we obtain the orthogonal polynomials P mi , m2 with respect to the constant measure \J/ (the 
polynomials are correspond to exited states) which satisfy the recursion relations (5-51). 



6 CTQ W on underlying graphs of two- variable P-polynomial 
association schemes via spectral method 

CTQW was introduced by Farhi and Gutmann in Ref.[5]. Let 1 2 {V) denote the Hilbert space of 
C-valued square-summable functions on V (i.e., J2% |/J 2 < oo). With each a E V we associate 
a ket | a), then {|a), a E V} becomes a complete orthonormal basis of 1 2 {V). 

Let \4>{t)) be a time-dependent amplitude of the quantum process on graph Y. The wave 
evolution of the quantum walk is 

ihj t \<t>{t))=H\<t>{t)), (6-62) 

where assuming H — 1, and \<f> ) be the initial amplitude wave function of the particle, the 
solution is given by \(/>o(t)) = e~' lHt \(f) }. It is more natural to deal with the Laplacian of 
the graph defined by L = A — D as hamiltonian, where D is a diagonal matrix with entries 
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Djj = deg(o;j) (recall that deg(<x,) is degree of the vertex a,- defined by the number of edges 
incident to the vertex a,-). This is because we can view L as the generator matrix that describes 
an exponential distribution of waiting times at each vertex. But on rf-regular graphs, D = ^7, 
and since A and D commute, we get 

e -itH = e -it{A-\i) = e -it/d e -itA^ (6 _ 63) 

this introduces an irrelevant phase factor in the wave evolution. In this paper we consider 
L = A = A\. Therefore, we have 

|0o(O> = e- M U>- (6-64) 
One approach for investigation of CTQW on graphs is using the spectral distribution method. 
CTQW on underlying graphs of P-polynomial association schemes has been discussed exhaus- 
tively in [15] via spectral method. In the following we investigate CTQW on underlying graphs 
of two-variable P-polynomial association schemes constructed in section 4 using spectral dis- 
tribution method. 

6.1 CTQW on underlying graphs of two- variable P-polynomial as- 
sociation schemes derived from Z m x Z m 

In the graphs constructed from Z m x Z m) the adjacency matrix is written as A z + A 2 and 
so we assign polynomial z + z to adjacency matrix. Then, by using the relation (5-50), the 
expectation value of powers of adjacency matrix A over starting site |0 O o) can be written as 

(0oo|^ m |0oo) = / j(z + z) m u.(z,z)dzdz, m = 0,1,2,.... (6-65) 

In the case of underlying graphs of two- variable P-polynomial association schemes, the adja- 
cency matrices are two-variable polynomial functions of A z and A z , hence using (4-46) and 
(6-65), the matrix elements (4>ki\A m \ O o) can be written as 

(<f>ki\A m | OO ) = -L(0oo|^ m | 0oo) = -^(<f>oo\Pki(A z ,A,)A m | O o) 
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]= I I (z + z) m P kl (z,z)u.( y z,z)dzdz, m = 0,1,2,.... (6-66) 
aid jr jr 



'dkl JRJR 

One of our goals in this paper is the evaluation of amplitudes for CTQW on underlying 
graphs of two-variable P-polynomial association schemes constructed in section 4 via spectral 
distribution method. By using (6-66) we have 

Pki(t) = (<f>ki\e- iAt \<l> o) = (<j>ki\<fa>(t)) = I I e- l{z +~ z)t Pki{z^{z,z)dzdz, (6-67) 

,/flH JR JR 



/dkl JRJR 

where (<f>ki\<f>oo(t)) is the amplitude of observing the particle at level kl (stratum r^((0, 0))) 
at time t. One should notice that, as illustrated in section 5, the polynomials pki(z,z) are 
obtained from realification of generalized Gegenbauer polynomials P m , n defined by (5-51). 
The conservation of probability J2k=o I (4>ki\4>oo{t)) | 2 = 1 follows immediately from (6-67) by 
using the completeness relation of orthogonal polynomials P mn (z,z). Obviously evaluation of 
{(f>ki\<f>oo(t)) leads to the determination of the amplitudes at sites belonging to the stratum 

r«((o,o)). 

Spectral distribution ji associated with the generators is defined as 

Z >3E^- z Kl )8{z - z^), (6-68) 
m k,i 

where k, I G {0, 1, m — 1}. Now using (6-67) and spectral distribution (6-68), the probability 
amplitude of observing the walk at stratum r i3 -((0, 0)) at time t can be calculated as 

P^t) = e" 2i '( cos ^ fe / m+cosM / m+cos2 ^( A;+ ')/ m )p ij (^ )Z , z* Kl ). (6-69) 

m k,i 

In particular, the probability amplitude of observing the walk at starting site at time t is given 
by 

-p>0o(0 \ \ ' — 2it(cos 27rfe/m+cos27ri/m+cos27r(fc+i)/m) (6-70) 

m2 k,l 

Example: case m — 3. 

By using (4-35), we obtain z M G {0, 3, 3a;, 3lu 2 }. Then by (6-68), spectral distribution is 
calculated as 

^z, z) = \{8{z - 3)5(z - 3) + 65(z)5(z) + 5{z - 3uo)5{z - 3a; 2 ) + 5(z - 3u 2 )5(z - 3u)}. (6-71) 
9 
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Therefore, by using (4-34) and (6-67), probability amplitudes of observing the walk at starting 
site, stratum r 10 ((0, 0)) and r n ((0, 0)) are calculated as 

P 00 (t) = ^{e- m + e 3U + Q}, 

P w (t) = 2 -{e~™ - e™}, 

P u (t) = \{e- m + 2e m - 3}, (6-72) 
9 

respectively. 

At the limit of the large to, we obtain the root lattice A 2 (hexagonal lattice). In the 
following we investigate the CTQW on root lattice A 2 using spectral method. 

6.2 CTQW on hexagonal lattice 

In this subsection we give continuous measure in the limit of the large size of the underlying 
graphs of symmetric two-variable P-polynomial association schemes derived by A z + A 2 in 
section 4. 

In the limit of large to, the roots z H = cu k + cu l + u~( k+l ^ reduce to z^i = e lXl + e tX2 + e~ l( - Xl+X2 ^ 
with x\ = lim^m^oo 2nk/m and x 2 = liim !m _ >0O 27tZ/to and the spectral distribution given in 
(6-68), reduces to continuous constant measure u,(xi,x 2 ) = 1/4-7T 2 . 

Also, the measure \i can be given in terms of complex variables z and z as 

H( z , z) = r dx x dx 2 b{z - (e ixi + e iX2 + e - l(xi+x ^))5(z - (e~ ixi + e~ ix ' 2 + e i{xi+X2) )) = 
Jo Jo 

= 1 1 (6-73) 

4n 2 ^-z 2 z 2 + A(z 3 + z 3 ) - 18zz + 27 

Then, the probability amplitudes Pki(t) are given by 

Pki(t) = (^|e- iA *|0oo) = r^i^e-^* 1 -^ 382 ^ 1 ^ 2 ^*,^!,^), (6-74) 

Jo Jo 
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where, Am = Pki(xi,X2)- In particular, the probability amplitude of observing the walk at 
starting site at time t, is calculated as 

P 00 (t) = f 2 " f 2 " dx l dx 2 e- m{cosxi+COSX2+cos{xi+X2)) . (6-75) 
jo Jo 

6.2.1 asymptotic behavior of quantum walk on hexagonal lattice 

As regards argument of the end of section 4.1, we can not obtain an analytic expression for the 
amplitudes of the walk in the infinite case, i.e., the integral appearing in the (6-74) is difficult 
to evaluate, but we can approximate it for large time t by using the stationary phase method. 
Studying the large time behavior of quantum walk naturally leads us to consider the behavior 
of integrals of the form 

/oo poo 
/ dx x dx 2 g{x)e- lt ^ x) (6-76) 
-oo J —oo 

as t tends to infinity. There is a well-developed theory of the asymptotic expansion of inte- 
grals which allows us to determine, very precisely, the leading terms in the expansion of the 
integral in terms of simple functions of t (such as inverse powers of t). Our basic technique 
will be to evaluate this integral in some approximation. The approximation we shall use will 
be the semiclassical expansion, which amounts to the well-known stationary phase approxi- 
mation as applied to the path integral. In this approximation, one can evaluate I(t) in (6-76) 
asymptotically as follows 

dx x dx 2 g[x)e-^ {x) ~ ]T g{d)e- itm — {DetA)- l/ \ (6-77) 

a 

where, summation is over all stationary points a of function f(x) and A is Hessian matrix 
corresponding to f(x). For more details about this approximation method, the reader is 
refered to [24]. 

Now, by using (6-77) we can discuss the asymptotic behavior of the amplitude of observing 
the quantum walk at starting site at large time t, i.e., we deal with the integral (6-76), with 
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g(x) = 1 and 

f(x) = 2(cosx! + cosx 2 + cos(x 1 + x 2 )). (6-78) 
Therefore, for the asymptotic form of the amplitude Poo(t) we get 

I (t) ~ -(i- e -«+«/2 + 1 at + 1 e 3<*-W2). (6.79) 

The asymptotic behavior of probability amplitudes Pfcz(t) (A;/ 7^ 00) at large time t, can be 
evaluated similarly. 

In order to obtain the asymptotic behavior of quantum walk on finite hexagonal lattice at 
large time t, we can compare the finite amplitude Poo{t) (Eq.(6-70)) and the continuous prob- 
ability (6-75) at large time t. Therefore, we calculate numerically the difference of amplitudes 
of the walk on root lattice A 2 with ones on finite hexagonal lattice, for large time t 

vr(m,t) = \I (t) -ij 1 e - 2 ^ os ^) +cos ^) +cos ^))|. (6-80) 
771 k,l=0 

The result has been depicted in Fig. 3. The figure shows that, the difference n(m,t) is 
limited to zero for m larger than ~ 50 and t ~ 1000. Therefore, to study the behavior of 
asymptotic quantum walk on finite hexagonal lattice, we can use arithmetic, approximate 
it with root lattice A 2 , and by using the stationary phase method, study the behavior of 
asymptotic quantum walk. 

6.3 CTQW on finite honeycomb lattice via spectral method 

In the graph T s constructed in section 4 from two copies of finite hexagonal lattices, we have 

A = 0+ ® B l + a- ® B, (6-81) 
where, B = I + S^ 1 + S 2 . Therefore, the spectrum of B is calculated as 

z u = 1 + uj~ k + J , k,l E {0, 1, ...,m - 1}; uj = , (6-82) 
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and the eigenvalues of A are given by 
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\ kl = ±\z M \ = ±^3 + 2(cos(2vr£;/m) + cos(27r//m) + cos(2tt(A; + l)/m)). (6-83) 

Now we apply spectral method by assigning variables z\ and z 2 to B and B f , respectively and 
a new variable z for er + and o"_ commonly. Then, the variable assigned to adjacency matrix A 
will be \zi\z + \z 2 \(z — 1). Clearly we have z\ = z 2 and so \z x \ = \z 2 \. Therefore, we will have 

A = \z 1 \(2z - 1). (6-84) 

The spectral distribution associated with adjacency matrix A is given by 

H{z x , zi, z) = ^n(zi, z^iSiz) + 5{z - 1)), (6-85) 

where, 

H(z u zi) = -^J] 6(z! - Zki)S(zi - z* M ). (6-86) 
m Tt 

Then, the probability amplitude of observing the walk at stratum 1^(0; 00) at time t is calcu- 
lated as follows 

Pk{t) = 1 1 / e " (|2ll(2 ^ 1)) V fe (^i;^(^i,^i^)^i^i^; k = 1,2,3,4, (6-87) 

where, A k = Pk(zi] z). In particular, the probability amplitude of observing the walk at starting 
site at time t is calculated as 

tee 1 m_1 / 

P (t)= / / / e-^ ZlK2z - 1))t fj,{z 1 ,z 1 ;z)dz 1 dz 1 dz = — ^ cos(y/3 + 2{cos2irk/m + cos2irl/m))t. 
J J J m k,l=0 

(6-88) 

For the sake of clarity, in the following we give details for the case m — 3. 
Example: case m = 3 
From the relations (4-39), we have 

A = l, A 1 =p 1 (z 1 ;z) = \z 1 \(2z-1), A 2 = p 2 ( Zl ; z) = (\ Zl \(2z - l)) 2 - 3, 
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A 3 = p 3 ( Zl ; z) = ±\ Zl \(2z-l)[(\ Zl \(2z-l)) 2 -5], A, = p 4 ( Zl ; z) = | [( | ^ | (2^- 1 ) ) 4 +3( | ^ | (2^- 1 ) ) 2 ] . 

(6-89) 

Using (6-82) and (6-83), the spectral distribution is obtained as 

= Mzi -3)<J(ii -3) + 25(z 1 -(2+uj))5(z 1 -(2+uj 2 )) + 25(z 1 -(2+uj 2 ))5(z 1 -(2+uj)) 
9 

+25(z 1 )5(z 1 ) + 6( Zl - (2u + l))S{z! - {2uo 2 + 1)) + S( Zl - (2uj 2 + l))<?(ii - (2cu + 1))}. (6-90) 

Using (6-88), the probability amplitude of observing the walk at starting site at time t is 
calculated as 

P (t) = ]- (cos 3t + 6 cos V3t + 2). (6-91) 
9 

Other probability amplitudes can be calculated using (6-87) and (6-89). 
6.3.1 CTQW on honeycomb lattice 

In the limit of the large m, the eigenvalues zm reduce to Zki = 1 + e~ tXl + e lX2 where x\ = 
lim^m^oo ^ and x 2 = liim^-^ ™. Therefore, the continuous spectral distribution is 

r2n p2tt 

p( Zl , Zi) = / / dx 1 dx 2 S(z 1 - (1 + e- ixi + e iiea ))<J(ii - (1 + e ixi + e"" 2 )) = 
Jo Jo 

/ 1 (6-92) 

4ttY3 - *W - (z\ + Z! 2 ) + 2z 1 z 1 (z 1 + z x ) - 2(zi + z x ) 

Therefore, in the limit of the large m, the probability amplitude of observing the walk at 
level k is given by 

P k (t)= (** rfxirfa;2(eV 3 + 2 ( COS2; i+ cos:E 2+cos(x 1 +x 2 ))t p ^ ;ri) ^ ^ +e ^z+2(co SXl +cosx 2+ co S (x 1 +x 2 ))t pk ^ X2 . l 
Jo Jo 

(6-93) 

In particular, for the probability amplitude Po(t) we have 

r p p27T r27T i 

Po(t) — J J cos(\zi\t)fj,(zii z~\)dz\dz~\ — J J cos y 3 + 2(cosx 1 + cosx 2 + cos(a; 1 + x 2 ))dxidx 2 

(6-94) 
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6.3.2 asymptotic behavior 

Similar to the finite hexagonal lattice, we investigate the asymptotic behavior of quantum walk 
at large time t using stationary phase approximation method. In this case we deal with the 
integral 

J Q (i^ = f [ e ~ i (y/3+2(cosx 1 +cosx 2 +cos(x 1 +x 2 )))t ( ^ x ^ x ^^_ f f g +i(Y / 3+2(cos xi+cos x 2 +cos(xi+x 2 )))t 

JO Jo Jo Jo 

(6-95) 

as probability amplitude Po(t), so, by using the stationary phase method we approximate the 
integrals for large time t. Here, we have g(x) = 1 and f(xi, x 2 ) = ^3 + 2 (cos a; i + cosx 2 + cos(:ri + x 2 )). 
Then the asymptotic form of the probability amplitude Po(t) is calculated as 

J„(t) ~ -(2VSsin3t + 2cost). (6-96) 

Now, we compare the finite probability amplitude (6-88) and the continuous probability (6-96) 
at large time t. Therefore, we calculate numerically the difference between amplitude of walk 
on honeycomb lattice and finite one, for large time t 

1 m_1 / 

n(m,t) = \I (t) C0S (V 3 + 2(cos2?rA; /m + cos2irl/m))t\. (6-97) 

171 k,l=0 

The result has been depicted in Fig. 4. The figure shows that, the difference n(m,t) is limited 
to zero for m ~ 60 and t ~ 700. Therefore, to study the behavior of asymptotic quantum 
walk on finite honeycomb lattice, we can approximate it with infinite one, and by using the 
stationary phase method, study the behavior of asymptotic quantum walk. 



7 Generalization to Z® n 

Similar to the case n = 2, we choose generating set P(io...o) = {(10...0), (0...01), (m — 
l...m — 1)} for Z® n . Then, the orbits of Weyl group S n+ i which are indexed by n-tuple 
m = (mi, ...,m n ), are given by 

P m = 0((mi, -m 2 , -m 2 - m 3 , -m 2 - m 3 - ... - m n )). (7-98) 
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By using (4-17) one can obtain a translation invariant partition R for (Z®™) 2 . In the regular 
representation, the adjacency matrix of underlying graph is written as 

A Z1 = S 1 + ... + S n + (SL.-Sn)- 1 , (7-99) 



with Si = I <S> ■■■ <S> J^®I <S> ... <S> I- Clearly, the relations R m define an abelian association 

i 

scheme (not necessarily symmetric) on Z^ 1 . Moreover, the corresponding Bose-Mesner algebra 
is generated by 

n+1 

A z k = J2 S h S i2 ...S ik , fc = l,...,ra, (7-100) 

i 1 <i 2 <...<i k 

where, S n+ i = (Si-.-Sn) -1 . In the regular representation of the group, for the corresponding 
adjacency matrices we have 

A m = J2 9 ■ (7-101) 

geP m 

From (7-100) and (7-101), it follows that the adjacency matrices satisfy the following re- 
cursion relations 

n+1 

A Zk A m = ^ A m+Vii+ ___ +Vi ^ , k — l,...,n, (7-102) 

n<i2<—<ik 

where, Vi, for i — 1, ...,n + 1, are n-dimensional vectors whose components are 

= , / = l,2,...,n. (7-103) 

In particular, A m = p m (A Zl , ...,A Zn ), where p m is a polynomial of degree mi + ... + m n with 
real coefficients. We refer to these types of association schemes as n-variable P-polynomial 
association schemes. 

Now, we symmetrize the graph to obtain an undirected underlying graph as in the case of 
n = 2. Then, we will have for the adjacency matrix 

A = A Zl +A Zn , (7-104) 

where, A Zl is given by (7-99) and A Zn = (AjJ*. The spectrum of A Zk (indexed by n-tuple 
I = (h, ...,/„)), is given by 

n+1 

z\ k) = Yl h G {0,l,...,m-l}, (7-105) 

i 1 <i 2 <...<i k 
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where, to = exp 2 ™/" 1 and l n+ \ = —(h + ... + /„). From (7-105), on can see that the spectrum of 
A Zn is complex conjugate of the spectrum of A Z1 . Therefore, the eigenvalues of A are given by 

^h,...,i n = 2(cos27r/i/m + ... + cos2irl n /m + cos27r(/i + ... + l n )/m). (7-106) 

The spectral distribution associated with the generators is given by 

*o = E - - - 4 n U (mot) 

where /, G {0, 1, m — 1}, for % = 1, n. Now using (6-67) and spectral distribution (7-107), 
the amplitude of observing the walk at level i = ...,i n ) can be calculated as 

p.(f) — \ -2it(cos27rii/m+...+cos27ri„/m+cos27r(ii+— +l n )/m) J (1) («) \ ('7-108") 

III i j 

where, the polynomials pi(zi, z n ) are given by the following recursion relations 

n+l 

z k P m = E Pm+v il+ ...+v ik ,k = l,...,n, (7-109) 

ji<j 2 <---<«fc 

where, Vi for 2 = 1, ...,n + 1 are defined by (7-103). In particular, the probability amplitude of 
observing the walk at starting site at time t is given by 

Pq(J^ ■= — !— g~ 2jt(cos27r/i/m+...+cos27ri n /m+cos27r(/i+...+/ n )/m) (7-110) 

In the limit of large m, the eigenvalues z\^..\ n reduce to z^.x n = e lXl + ... + e tXn -\- e -< x ^+-+ x ^) 
with Xi = \imi itm ^ 00 2nli/m, and the spectral distribution reduces to continuous constant 
measure fj,(xi, ..,x n ) = l/(27r)™. In fact, in the limit of large m, the study of CTQW on finite 
symmetric graph constructed from Z® n as above, is equivalent to the study of walk on the 
root lattice A n , where the continuous form of probability amplitude Pi(t) is given by 

Pi(t) = r ... / 2 %-^(co S , 1+ ... + co Sa;rl+ cos(, 1+ ... + a ;rl )) pi(xi) ... >a . n ) da . 1 ... da . n> (Mil) 

{Any 1 Jo Jo 

where, Ai = Pi(x±, x n ). In particular, the probability amplitude of observing the walk at 
starting site at time t is given by 

P (t) = — 1— f 2W ... / 2 " e - 2i ' (cos:ci+ - +cos:E " +cos{a;i+ - +:E " )) rfxi...rfa; n , (7-112) 
(27r) n Jo Jo 
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where, the integral in (7-112) can be approximated by employing stationary phase method at 
large time t. 

8 conclusion 

Using the spectral distribution method, we investigated CTQW on root lattice A n and honey- 
comb one, by constructing two types of of association schemes and approximating the infinite 
lattices with finite underlying graphs of constructed association schemes, for large sizes of the 
graphs and large times. Although we focused specifically on root lattice A n and honeycomb 
one, the underlying goal was to develop general ideas that might then be applied to other 
infinite lattices such as root lattices B n , C n and etc. also quasicrystals with certain sym- 
metries. Apart from physical results, we succeeded to obtain some interesting mathematical 
results such as a generalization to the notion of P-polynomial association scheme, where we 
expect that, the n-variable P-polynomial association schemes posses the analogous properties 
of P-polynomial association scheme and can be applied in coding theory in order to construc- 
tion of new codes. We hope that, studying other infinite lattices leads us to other interesting 
mathematical objects, perhaps new types of association schemes. 

Appendix : Suppersymmetric structure of T s 

From the block form of A in (4-37), we can see that the constructed underlying graph of 
association scheme from two copies of finite hexagonal lattice, has supersymmetric structure. 
Following Ref.[30], we introduce the model of supersymmetric algebra as follows 
We define two operators Q + and Q- as 



O O 



( 



O B f 



Q+ 



Q- = 



(A-i) 



V 



B O 



) 



\ 



o o 



J 



(A-ii) 
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Then we define two hermitean charges Qi, Q2 and Hamiltonian H as follows 

Qi=Q+ + Q~, Q 2 = -i(Q+-Q-), H = Q\ = Q\. 
With the above definitions, we get 

Ql = Q l = 0j H = {Q + ,Q^}, [H,Q±]=0, [H,Q lt2 ]=0, 

{Qi,Q 2 } = 0, =s> {Q i ,Q j } = 2H5 ij . 

Therefore, Q + , Q- and H generate a closed supersymmetric algebra. 

For the association scheme derived from two copies of finite hexagonal lattice in section 4, 
we make the following correspondence 



(A-iii) 
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(A-iv) 



and therefore, 



H = A 2 
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(A-v) 



In other words, the adjacency A is our original Dirac operator. Now it can be checked that all 
the above (anti) commutation relations are fulfilled by our representation in the form of graph 
operators. In our special graph, B and B l commute with each other, so the Hamiltonian is of 
the form H — I <S> B l B. Therefore, the spectrum of H is at least twofold degenerate, i.e., 



H(f,g) t = E(f,g) t 



B l Bf = Ef, B t Bg = Eg. 



(A-vi) 



Hence, (/, 0)* and (0, g) 1 are eigenvectors of H to the same eigenvalue. 

As H = Q\ = Q\ and {Qi,^} = 0, certain combinations of the above eigenvectors yield 
common eigenvectors of the pairs H, Qi. 

From the fact that [B, B l ] = 0, we know that B and B l have common eigenvectors. If 



Bf=Xf, B t f=Xf, 



(A-vii) 
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then 

BB f f = XX' f. (A-viii) 
As the spectrum of B is complex conjugate of the spectrum of B l , we have 

BB l f= \X\ 2 f. (A-ix) 

In other words, the eigenvector of B with eigenvalue A, is eigenvector of Hamiltonian H, with 
eigenvalue |A| 2 and degeneracy at least 2. 
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Figure Captions Figure-1: Shows root system corresponding to A 2 . 

Figure-2: Shows finite honeycomb lattice with generators cl, cS'f 1 and cS 2 . 

Figure-3: Shows n(m : t) for root lattice A 2 as a function of m (where, m 2 is the 
number of vertices of finite hexagonal lattice) at t ~ 1000, where the difference is 
almost negligible for m > 50. 

Figure-4: Shows 7i(m,t) for honeycomb lattice as a function of m (where, 2m 2 is 
the number of vertices of finite honeycomb lattice) at t ~ 700, where the difference 
is almost negligible for n > 60. 
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